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Abstract 

We show that within classical statistical mechanics it is possible to naturally derive power law 
distributions which are of Tsallis type. The only assumption is that microcanonical distributions 
have to be separable from of the total system energy, which is reasonable for any sensible 
measurement. We demonstrate that all separable distributions are parametrized by a separation 
constant Q which is one to one related to the q-parameter in Tsallis distributions. The power-laws 
obtained are formally equivalent to those obtained by maximizing Tsallis entropy under q 
constraints. We further ask why nature fixes the separation constant Q to 1 in so many cases 
leading to standard thermodynamics. We answer this with an explicit example where it is possible 
to relate Q to sytem size and interaction parameters, characterizing the physical system. We 
argue that these results might be helpful to explain the ubiquity of Tsallis distributions in nature. 
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INTRODUCTION 



There has been a tremendous interest in a generahzed definition of entropy, recently 
introduced by TsaUis 

where p is the the normahzed energy density and dT indicates phase space integration. The 
reason why this modification has attracted so much interest is partly because of the possibil- 
ity to derive power-law distributions in the canonical ensemble within the maximum entropy 
principle. This modification of entropy and its resulting formalism, which is sometimes re- 
ferred to as non-extensive thermostatistics, has triggered far more than a thousand works 
in the past few years Q]- However, despite its phantastic descriptive success of power-laws 
in physical, chemical, biological, and social systems, it has not yet been possible to derive 
this form of entropy from thermodynamic or statistical principles. Within the formalism 
suggested by Tsallis it is necessary to define expectation values of quantities depending on 
the energy spectrum of the system not in the standard way but with so-called o-expectations 

nn ri 

recover the Legendre structure of thermodynamics p|. The physical 
interpretation of these g-expectations is under heavy debate. 

There have been several papers recently with the aim to derive canonical power distribu- 
tions from first principles, see e.g. 0,Q,0- In convincing work it has been beautifully 
noted that the expression ^ ^-^ j = g — 1 gives a physical meaning to q and that power laws 
in the canonical ensemble can be derived on a Hamiltonian basis. 

In this work we adopt a somewhat different strategy and derive - using a mathematical 
theorem - power-law distributions for the canonical ensemble directly, just by the use of the 
variational principle, and a separation Ansatz, without touching the standard definition of 
entropy. 



POWER-LAWS IN THE CANONICAL ENSEMBLE 

We begin by noting that any thermodynamic system which can be measured in equi- 
librium must be separable, i.e., that the thermodynamic quantities of the measured system 
should not explicitely depend on the energy of the total system E. In the following we 
consider a sample (observed system) in contact with a reservoir. The energy of the sample 
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is El, the energy of the reservoir is E2, such that the total (isolated) system has a constant 
total energy E = Ei + E2. The number of microstates in the sample is u!i{Ei) and uj2{E2) in 
the reservoir. The energy of the sample fluctuates around its equilibrium (extremal) value 
denoted by E^. The Hamiltonians describing the sample and the reservoir are Hi and 
respectively. Thermal contact of the two systems means H = Hi + H2 and the partition 
function Z{E) is the convolution of the two microcanonical densities 

E 

Z{E) = j dEiUi{Ei)u2{E - El) , (2) 


with 

uj,{Ei) = J dV, 6{H, - Ei) . (3) 

Following the usual line of thought to pass from the microcanonical to the canonical de- 
scription, represented by p is given (up to a constant multiplicative factor) by 

p{Ei) = uji{Ei)uj2{E - Ei)Z-\E) . (4) 

Note, that this description is dictated by the equations of motion. Assuming the existence 
of a unique extremal configuration at some Ei = E^. defined by 5p = 0, leads to the well 
known condition 

UJi^Er=E. ^^\e,=E-E,- ^ ' 

which defines the temperature T of the system. The usual definition of entropy Si = k ln{uJi) 
implies that the extremal configuration is found where S = Si + S2 is extremal with its 
associated temperature as defined above. 

Under which circumstances can one factorize the dependence of p on the total energy 
E7 We are hence looking for classes of microcanonical distributions that allow for such a 
separation of E into a multiplicative factor. A standard way to motivate the appearance 
of the Boltzmann term in the canonical ensemble can be seen as a consequence of this 
i?-separation 

iU2iE-Ei) = exp{\n{uj2{E - El))) 

^ exp{\niu;2{E))~^\n{LU2)Ei) (6) 

^ U2{E)exp{-PEi) . 
It is worthwhile to note that the approximation in Eq. Q is exact for uj2{E — Ei) being 
an exponential in E. Up to this point we have summarized textbook knowledge. It is one 
purpose of this work to emphasize that (jH)) is not the most general way of separation. 



To find the most general separation, we generalize the log function in Eq. © to some real 
function /, being strictly monotonous and twice differentiable. Monotonicity is needed for a 
well defined inverse f-K The idea is to write El) = f-^ o f ouj{{E - E^) - {Ei~ E^)) 

and to expand f o uj around E — E^,. Suppose energy E is separable from the system, then 
there exist two functions g and h such that 

u;{E-E,) = g{uj{E-E,)) h{x) , (7) 

with X := l3{Ei — -E*); to simplify notation we write u := uj{E — E^) in the following. We 
now use / to find the unknown functions g and h by expanding f o u to first order 

f{u;{E-E,))=f{g{u;)h{x))^f{u;)-u;xf{u;) , (8) 

which is justified for small x, i.e., the system being near equilibrium. The most general 
solution to this separation Ansatz is given by the family of equations (/, g, h)Q, parametrized 
by a separation constant Q, and C and C2 being real constants 

/H = c^i-« + C2 

9{uj) = u; (9) 
h{x) = [1-{1-Q)x]^ 

To see this, first g is found by setting x = and ho = h{0), so that Eq. (jH)) yields 
f{g{Lj)ho) = f{oj), which means g{oj) = Without loss of generality set ho = 1 and 
arrive at f{u!h{x)) = f{uj) — ujxf'{uj). Form partial derivatives of this expression with 
respect to x and a), and eliminate the f'{ujh) term from the two resulting equations 

f{uh)h' = -f'{u) ^^^^ 

f{uh)h = {i-x)r-ujxf" 

to arrive at the separation equation 

where Q is the separation constant. The differential equation 1 ~ + l) = <5 is straight 

1 

forwardly solved to give h[x) = [l — (1 — Q)x'j using h{0) = 1 to fix the integration 
constant. The equation —Ujj^ = Q means, f{iu) = Cij^uj^^'^ + C2, with Ci and C2 
integration constants. / is strictly monotonous except for Q = 1, where it is constant. It is 
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straight forward to test that Eqs. © solve Eq. (jH)). The term of interest in the canonical 
distribution can now be written as a generalized Boltzmann factor 

uj2{E-E{)=uj2{E-E,)[l-{l-Q)P{Ei~E,)\^ . (12) 

The separation constant is not specified at this level. As we will see below the choice of 
a particular physical system will determine Q. The usual Boltzmann factor Eq. © is 
recovered as a special case in the limit Q ^ 1. Note, that if uj2 is of the form uj2 oc E^^^^'^, 
then Eq. (fT^ holds exactly and not only to the first order approximation in Eq. (jHl). The 
best way to prove this is to write /3 = ^\e-e, = (^i_q-)i^e-e*) compute straight 

forwardly 

u;2iE - E,) = {E- E,)T^ (l - fff ) ^ 
= UJ2{E-E,)[l-{1-Q)I3{E^~E,)\^ , 

see also llll . As we will see below, uj2 covers classical (homogenous) Hamiltonians 

with pair-potentials. 

Having derived the principal form {f,g, h)Q to first order it is easy to see that with the 
same family of functions we can expand 002 to all orders. Repeated differentiation of p 
at its extremum leads to a hierarchy of equations relating properties of the u densities to 
properties of p at its extremum. With the definitions 



n 



[n] being the n th derivative, and 2 = 1, 2 indicating system 1 and 2 {Eu = E^,, E2* = E—E^). 
We successively construct the whole hierarchy to find 

r2 = 4>l + 4>l-2 

rs = (0^-30i)-(0i-3(/.i) 

U = (01 - m + 30i(r2 + 2 - 01)) (15) 

+ (02 _402 + 302(^^ + 2-02)) 



Using these equations we can re-express the 0^ terms in terms of (f>\ and r„, so that all we 
know about system 2 is encoded in local properties of p at the equilibrium. To fourth order 
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the general expansion reads 

Lu{E - E^) =iuiE- E,){l + il-Q)[- 
(16) 



X 
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+ l,(cj)l-Q)x'-l,{-3Q<Pl + cPl + {Q + l)Q)x 



- {Q + 2){Q + l)Q)x' + ...]}^ , 
which in the hmit Q 1 is 

uj{E - E^)=uj{E~E,)exp{-x + ^{(j)l-l)x^ 

- M-30i + 01 + 2)x3 + 1(01 -401 -3(01)2 (17) 
+ 1202-6)x^ + ...} 

Note, that for uj2 oc exp{i3E) all 0^ = 1, and all higher order terms vanish and the standard 
Boltzmann result becomes exact. This concludes the main finding of the paper. 



THE PHYSICAL MEANING OF THE SEPARATION CONSTANT Q 

Why does nature fix Q — > 1 in so many cases, i.e. why is standard thermodynamics 
the most predominantly realized situation? In the following we will demonstrate with the 
help of an example how the separation constant Q appearing in Eq. ()12j) can be related to 
system size and interaction parameters of a real physical system. This explains the ubiquity 
of Q = 1. Examples of this kind have been given in a somewhat different context before 



An example 

Let us specify the following particle Hamiltonian for pair-potentials, governing the 
sample in D space dimensions. We use n = DN . 

N 2 ^ 

^(^'^^) = E|^ + Ei^^-^^r ■ (18) 

i i<j 

The energy density is given by the phase space integral 

,{E) = jd-pd-x 5 (e. It + Ek, i^^ - ^.r - e) 

= J^'^dE,Jd-pd-x6{E.£-E,) (19) 



UJ 



We compute the kinetic term 



V 2m 



n-l 

and the potential contribution 



(20) 



X 



k\2 



J (T'x S{E 



E't jd^'y 5{E 



Sj=2 Y^i=i Tl,k=i [vi ~ Vj) 



(21) 



E a 



-1 



const. 



where we used the substitution yi 
finally get for Eq. (fTTIjl 



X 



i/E^ and the fact / dx5{Xx) = J dx X ^6{x). We 



UJ 



[E) oc / dEx El [E-Ei)" oc E 
Jo 



(22) 



This allows us now to compare exponents (and coefficients) in Eq. ()22|1 and Eq. (fT^ to 
arrive at the relation 

' + 1 , (23) 



1-Q 2a 

which fixes the separation constant. This equation establishes the connection between the 
interaction term in the Hamiltonian a, the dimensionality of the phase space n = DN, and 
the separation constant Q. From Eq. ()2Hj) it is immediately clear that for large systems 
the separation constant is always Q 1, i.e. the classical Boltzmann term (jHl) is recovered. 
For small systems, with a fixed number of particles Q depends on the interaction between 



the particles. For an ideal gas a 



-00 the separation constant is Q = 1^^. Nontrivial 



Q ^ 1 should be expected for strongly interacting and/or small systems, i.e., |a;|/|a; + 2| ~ n 
and systems with —2 < a < 0, where the limit n —>■ 00 implies n{a + 2)/(2a) — 1 ^ —00 
so that BG is not obtained. Let us assume that for a system Q 7^ 1 is due to its finite 
size. Here the standard argument in BG thermodynamics of extremely sharp peaks of 
the distribution at equilibrium is not necessarily valid, and E^: and the expected energy 
U = {H) may significantly differ. Assuming uJi{Ei) oc E^' and the validity of the variational 



principle Sp = { 



j^dE,Et{E-E{) 




With this we can write 



(24) 



with (x) 



da;a;(ni+x)"i(n2-x)"2 



. Rewriting leads to 



(25) 



U 



n2-{x)\2 
ni + {x)i2 



1 + 



One can therefore apply the same equilibrium argumentation also to the (measurable) ex- 
pectation values by substituting ni = ni + (x)i2, n2 = '"■2 — (a^)i2 and p oc E^^{E — EiY^ so 
that 5p = with 5E = 0. The Q of the observable canonical ensemble is then 



for the example of a pair potential Hamiltonian for the reservoir. 
CONCLUSION 

Equation p2j] is without doubt formally equivalent to the Tsallis q form of the Boltzmann 
term. Our result (fT^ appears as a consequence of the separability of the total system 
energy E from the microcanonical density of the sample system with the parameter Q being 
nothing but a separation constant not further specified at this stage. Q is fixed and a 
physical meaning is obtained as soon as a particular Hamiltonian is specified. Equation (fT^ 
is not exactly what Tsallis gets by extremization of Sq, with his constraints for the canonical 

^-^ q 

ensemble: — 1; — y^^'?' = canst. In our derivation we get E^, and the 

classical (3 (real temperatures), where Tsallis gets the terms Uq and P/^iVh respectively. 
Our general Boltzmann factor is obviously not obtained as a result of a maximization of 

We have computed higher order terms and found that if the canonical distribution is a 
power, all higher terms in the expansion vanish and our result ()12|1 holds exactly and not only 
as a first order approximation as we started out with in (jHl). This might be an interesting 
finding since many physically interesting microcanonic densities behave like powers in Ei, 
at least for finite size systems. 



1 



(26) 



1-g 



~- = n2 
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An important aspect in our work is that all our arguments are strictly based on Hamilto- 
nians and on the variational principle. At no point we are forced to take the thermodynamic 
limit. For that case we had some discussion if the equilibrium point is different from the 
expectation value. In case that mean and equilibrium do not coincide, our result would of 
course not be observable, since measurements take place at the expectations whereas we 
have developed our results around the extremal configuration. We can however show that 
for microcanonical distributions characterized by u!i{Ei) oc E'^^ and uj2{E2) oc E'^'^, mean 
and extremum coincide for ki = k,2- In the thermodynamic limit the difference is of course 
completely irrelevant. 

The purpose of the given example is not to relate the first part of the paper to any non- 
extensive system, but only to demonstrate the possibility to relate the separation constant 
to physical parameters. However, we think that it is in principle possible to use the central 

jart of the paper to think about the occurrence of phenomena in situations described in e.g. 

1^ . Here within the long meta-stable regions we have the situation of having an 'almost- 
equilibrium' which might be characterized by a non-trivial Q before the t oo limit is 
taken. 

We thank H. Grosse for his interest and C. Tsallis for useful comments. 
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